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ABSTRACT 



Equilibrium molecular dynamics simulation and the Green-Kubo formalism 
were used to calculate self-diffusion coefficient, shear viscosity, and thermal 
conductivity for 38 different dipolar two-center Lennard- Jones fluids along the 
bubble line and in the homogeneous liquid. It was systematically investigated 
how anisotropy, i.e. elongation, and dipole momentum influence the transport 
properties. The reduced elongation L* was varied from to 1 and the reduced 
squared dipole momentum was varied depending on the elongation as follows: 
for spherical fluids (L*=0) from fi*^^0 to 20, for L*=0.2 from i^*^=0 to 16, 
and for L*=0.4 and above from iJ,*^—0 to 12. This represents the entire range 
in which parameters for real fluids are expected. The statistical uncertainty of 
the reported data varies with transport property, for self-diffusion coefficient 
data the error bars are typically lower than 3 %, for shear viscosity and thermal 
conductivity they are about 8 and 12 %, respectively. 

KEYWORDS: Green-Kubo; molecular dynamics; dipole; self-diffusion; shear 
viscosity; thermal conductivity. 
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1 INTRODUCTION 



In an accompanying work [IJ we report on a comprehensive investigation 
on self-diffusion coefficient, shear viscosity, and thermal conductivity of two- 
center Lennard- Jones plus point quadrupole (2CLJQ) fluids by equilibrium 
molecular dynamics and the Green-Kubo formalism. In the present work, 
that investigation is extended to two-center Lennard- Jones plus point dipole 
(2CL JD) fluids using the same methodology. These two model types are aimed 
at molecules with different type of polarity, i.e. 2CL JQ models can well be used 
for quadrupolar fluids such as CO2 or Chlorine [2], whereas 2CLJD models 
are for dipolar fluids such as CO or partly fluorinated alkanes [3]. 

Static thermodynamic properties of the 2CLJD fluids have been extensively 
investigated in the past. There are many results for vapor-liquid equilibria 
[ill5]|6irfll8] . excess properties [9], or surface tension p^pT] . The 2CLJD poten- 
tial has been successfully applied to modeling real fluids, yielding good results 
for vapor-liquid equilibria [T2|13|ll4f3] . Joule- Thomson inversion curves |15j . 
virial coefficients [TB], and recently also for shear viscosity and thermal con- 
ductivity of refrigerants [T7] . On the other hand, transport properties for this 
useful model class have not been explored systematically. An exception is the 
work of Lee and Cummings [T8] on shear viscosity using the Stockmayer poten- 
tial [in], which is a subgroup (L*=0) of the 2CLJD model fluids investigated 
here. The present work aims to fill this gap. 

This investigation, as the accompanying one on 2CLJQ fluids, is carried out 
for state points along the bubble line and in the homogeneous liquid. It covers 
38 model fluids, specified each by a certain combination of elongation and 
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dipole momentum. Thermodynamic conditions were selected on the basis of 
previous work of Stoll et ah [8)120] on vapor-hquid equihbria of 2CLJD fluids. 
As the present work is analogous to the accompanying one on 2CLJQ fluids 
[l], redundant information like the Green-Kubo equations is omitted here. 

2 MOLECULAR MODEL 

The intermolecular interactions are represented by the two-center Lennard- 
Jones plus point dipole (2CLJD) potential. The 2CLJD potential is pairwise 
additive and consists out of two identical Lennard- Jones sites a distance L 
apart (2CLJ) plus a point dipole of momentum /i placed in the geometric 
center of the molecule, oriented along the molecular axis connecting the two 
Lennard- Jones (LJ) sites. It is analogous to the quadrupolar 2CLJQ potential 
investigated in the accompanying work [T], only differing in the polar interac- 
tion that is a point dipole instead of a point quadrupole. The contribution of 
a point dipole is given by [21] 



wherein Yij is the center-center distance vector of two molecules i and j. 6i is 
the angle between the axis of the molecule i and the center-center connection 
line and (pij is the azimuthal angle between the axis of molecules i and j. 
Finally, eo is the electric constant 8.854187817-10"^2 ^2^^ 

Analogously to specific 2CLJQ models, a 2CLJD model for a real substance 
is fully determined by five parameters: a, e, L, [8] and the molecular mass 
m. In the reduced form also for this model class only two molecular param- 
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eters remain, i.e. reduced elongation L* = Lja and reduced squared dipole 
momentum /x*^ = /i^/(47reoe(j^). Henceforth, "squared" will be omitted in the 
text for brevity. 

Furthermore, all results of this study were obtained and are presented in the 
reduced form; the definitions are given in the accompanying work [1]. For the 
sake of brevity, " reduced" will be omitted in the following. 

3 INVESTIGATED MODELS AND STATES 

In the present work, 38 different model fluids were studied, where each fluid is 
fully determined by one combination of elongation L* and dipole momentum 
/x*^. Simulations at 12 liquid state points were carried out for each model fluid. 

The studied model fluids have elongations that vary from L*=0, i.e. spheri- 
cal molecules, to L*=l, i.e. strongly elongated dumbbell- shaped molecules, in 
seven steps. For Stockmayer fluids, seven dipole momenta from /i*^=0 to 20 
were considered. For fluids with L*=0.2, six dipole momenta from /i*^=0 to 16 
and for larger elongations five dipole momenta from /i*^=0 to 12 were taken 
into account. The upper limit of 12 is sufficient to describe strongly dipolar real 
fluids, e.g. CClFs with /i = 1.8261 D (/x*^ = 3.4932), CClsF with /x = 2.7009 
D (/x*2 = 3.6250), or CH2F-CF3 with /i = 3.0214 D (/x*^ = 8.0004). Very high 
dipole momenta (/i*^=16 and 20) were only considered for fluids with a small 
elongation, because they are realistic only for such molecules [3]. 

Due to consistency, spherical fluids (L*=0) were treated as two superimposed 
Lennard- Jones sites. To compare with one-center Lennard- Jones results, the 
temperature has to be divided by 4 as well as the dipole momentum. Also the 
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present values of self- diffusion coefficient D*, shear viscosity 77*, and thermal 
conductivity A* have to be divided by 2. 

To facilitate a meaningful comparison of the transport properties, another 
reduced form of the temperature Tr=T*/T* and the density Pb=P*/Pc w^ls 
used, where T* is the critical temperature and p* the critical density of the 
individual 2CLJD fluid; values for T* and p* were taken from [8]. 

For each fluid, three state points along the bubble line were studied from 
Tr=0.6 to 0.9 with temperature increments of ATr=0.15. In addition, an- 
other nine state points were studied in the homogeneous liquid on three iso- 
chores, starting from these three bubble points with temperature increments 
of ATr=0.15, cf. Fig. [H This grid is less dense than the one used in the ac- 
companying work on quadrupolar fluids [1] , however, as more model fluids are 
regarded here, both present roughly the same number of data points. 

For quadrupolar fluids critical and bubble densities in terms of p* increase 
monotonously with increasing quadrupole momentum [8], but dipolar fluids 
behave different [3]. Fig. [2] illustrates the variation of the bubble density, for 
three reduced temperatures as a function of the dipole momentum for Stock- 
mayer fluids (L*=0). As can be seen, the critical density, i.e. bubble density at 
Tr=1, increases up to /i*^ ~ 3, but decreases for higher dipole momenta. At 
Tr=0.9 and below they are larger than those of non-polar fluids in all cases, 
where for this temperature a significant maximum of the bubble density at 
intermediate dipole momenta is present. These tendencies are also found for 
elongated fluids. The influence of elongation on phase behavior does not differ 
from that of quadrupolar fluids [1], i.e. bubble densities decrease monotonously 
with increasing elongation [8]. 
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4 SIMULATION DETAILS 



As many of the considerations and technical details of the simulations are 
the same as those used in the accompanying work [Ij, we proceed describing 
only the differences. Here, the significant long range contributions due to the 
dipolar interaction were treated with the reaction field technique [23f2^ . Sim- 
ulations were carried out in the canonical ensemble using the modified equa- 



tions of motion with a time step of Aty e/m/cr = 0.001, and the Nose-Hoover 
thermostat P^PT] with a thermal inertial parameter of 10 in reduced units. 
Statistical uncertainties of the transport coefficients were estimated using the 
standard deviation of three independent simulations runs with 5 000 indepen- 
dent autocorrelation functions. Independence between correlation functions 
was achieved with a time span of 0.2 in reduced units between consecutive 
correlation functions. The convergency issues of self-diffusion coefficient [2H] 
and shear viscosity [29|30|I31|I32|33] are discussed in the accompanying work \T\ . 
The integrals usually converge within their statistical uncertainties at about 
1 in reduced units, nevertheless an integration time of At* =2 was used. 

5 RESULTS 

In this section, the simulation results for the transport coefficients are pre- 
sented. Numerical data for self-diffusion coefficient, shear viscosity, and ther- 
mal conductivity are given in Table [T] for spherical fluids (L*=0) with dipole 
momenta from fi*'^=0 to 20, for slightly elongated fluids (L*=0.2) with dipole 
momenta from /i*^=0 to 16, and for elongated fluids from L*=OA to 1 with 
dipole momenta ranging from /i*^=0 to 12. All data in Table [T] correspond to 
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state points along the bubble line for reduced temperatures of Tr =0.6, 0.75, 
and 0.9. The complete data set, with nine additional state points in the liquid 
region for each fluid, is available in [3l] and partially included in Figs. El El [HI 
l9l [TTl and [IS The effects of elongation, dipole momentum, temperature, and 
density are discussed in the following for each transport coefficient separately. 

The accuracy of the calculated transport properties decreases in the sequence 
self- diffusion coefficient, shear viscosity, thermal conductivity. The high accu- 
racy of the self-diffusion coefficient, with error bars lower than 3 %, is due to its 
individual nature [32] ■ Shear viscosity and thermal conductivity show larger 
uncertainties, that are around 8 and 12 %, respectively. In most simulations of 
the present work the autocorrelation functions of thermal conductivity decay 
faster than those for shear viscosity, but fluctuate more. 

Other factors that influence the accuracy of the reported data are elongation 
and dipole momentum. In particular, at low temperatures, for fluids with large 
anisotropy and strong dipole momentum, the transport coefficients show larger 
simulation uncertainties. 

In the following, the results are discussed for ten selected fluids, covering the 
whole range of the two molecular parameters, from spherical (L*=0) over elon- 
gated (L*=0.505) to strongly elongated (L*=1.0) fluids with varying dipole 
momentum of /i*^=0, 6, 12, and 20. A subset of six fluids is taken in some 
cases only due to graphical reasons. 
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5.1 Self- diffusion coefficient 



Figs.[3]and|l]illustrate the self-diffusion coefficient along the bubble line for ten 
selected fluids. The results can either be discussed in terms of reduced density 
Pji as in Fig. [HJor in terms of number density p* as in Fig. |H From Fig. [H] it can 
be seen that the regarded range of reduced density is similar for all fluids, but 
significant deviations from the principle of corresponding states are present 
also for the density. At constant Tr, it can be discerned that the self-diffusion 
coefficient always decreases with increasing elongation. The dipole, however, 
can either decrease or increase the self-diffusion coefficient. A better visibility 
of the data (which is even more needed for the less accurate properties shear 
viscosity and thermal conductivity) is obtained when plotted over number 
density in Fig. |H Therefore, this graphical representation is preferred in the 
following. 

As Fig. m shows, D* decreases with increasing number density along the bub- 
ble line. Self-diffusion coefficients of fully elongated fluids lie roughly along the 
same line as observed for quadrupolar fluids ^J. To analyze the effect of the 
dipole momentum, it is helpful to compare the bubble densities of Stockmayer 
fluids in Fig.[2l where the influence is most visible, with the corresponding self- 
diffusion coefficients in Fig. HI The self-diffusion coefficient shows a minimum 
for Stockmayer fluids with /i*^ f» 6 at all three temperatures. The bubble den- 
sities of Stockmayer fluids also show a peculiar behavior: maxima at Tr=0.9 
and 1 and a point of inflexion at Tr=0.6. As the self-diffusion coefficient of 
non-polar fluids decreases with increasing density, it can be concluded that 
the isolated effect of the dipole is to increase the self-diffusion coefficient. 
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Fig. O shows the dependence of D* on number density in the homogeneous 
hquid region at a constant reduced temperature of Tr=0.9. Note that the 
density range is the same as in Fig. HI Along this isotherm D* decreases slightly 
hyperbolic with increasing density, resembling the behavior of D* along bubble 
lines for a given elongation. Comparing D* along bubble lines with isothermal 
data for the same density variation, it is found that density is the dominant 
variable, being responsible for about 80 % of the variation in D*. 

Fig. [H] shows the dependence of the self- diffusion coefficient on reduced tem- 
perature for a subset of six selected fluids at different isochores. The isochores 
correspond to bubble densities at the reduced temperature Tr=0.6, cf. Fig. [H 
which have similar values in terms of pr. It can be seen that the self-diffusion 
coefficient decreases monotonously with increasing dipole momentum for elon- 
gated fluids, whereas for Stockmayer fluids, D* again exhibits a minimum at 
intermediate dipole momenta. Along an isochore, the self-diffusion coefficient 
increases linearly with increasing temperature. The gradients with respect to 
reduced temperature are almost constant for a given elongation but less pro- 
nounced for more elongated fluids. Such a behavior was also found in the 
accompanying work on quadrupoles [Ij. 

5.2 Shear viscosity 

Fig. [7] illustrates the shear viscosity along the bubble line for the ten selected 
fluids. In contrast to quadrupolar fluids, shear viscosity results for fluids with 
a given elongation, but different polar momentum, are not along a single line. 
At constant Tr, shear viscosity increases considerably with increasing dipole 
momentum. In accordance to quadrupolar fluids pLj, more elongated molecules 
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have generally a lower shear viscosity along the bubble line. The extremely 
high values of shear viscosity found for high dipole momenta (/i*^= 20) at low 
temperatures are remarkable. 

Fig. [8] illustrates the isolated effect of density on shear viscosity for the same 
ten fluids in the homogeneous liquid region at Tr=0.9. Comparing the varia- 
tion of ?7* along bubble lines and along isotherms in the same way as for D* , 
it is found that the density effect is responsible for about 80 % of the increase 
of 7]* along the bubble line. Exceptions are the shear viscosities of strongly 
dipolar Stockmayer fluids with /i*^=16 and 20. These are very sensitive to 
temperature changes, where the contribution of density is only about 10 % for 
/i*2=20. 

Fig. [9] shows the dependence of shear viscosity on reduced temperature for a 
subset of six fluids along different isochores with similar values in terms of pr. 
As expected, the shear viscosity decreases with increasing temperature. Non- 
polar fluids show little temperature dependence. On the other hand, strongly 
dipolar fluids, with an about twofold higher shear viscosity than non-polar 
fluids in the cold liquid, are more sensitive to temperature exhibiting larger 
gradients. 

5.3 Thermal conductivity 

Fig. [TU] illustrates the thermal conductivity along the bubble line. The thermal 
conductivity increases with increasing density and more elongated molecules 
have lower thermal conductivities, similar to quadrupolar fluids [1] . In contrast 
to quadrupolar fluids, the data for a constant elongation do not lie on a single 
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line here. It can best be seen for Stockmayer fluids that A* increases strongly 
with increasing dipole momentum at constant Tr. 

Fig. [TT] shows the density dependence of the same ten fluids in the homo- 
geneous liquid at T/j=0.9. As can be seen, the curves resemble those along 
bubble lines, cf. Fig. [TOl demostrating the dominant density effect. Similar 
results were found for quadrupolar fluids [1]. 

Fig, shows the weak temperature dependence of A* through isochoric data 
for a subset of six fluids. Taking the statistical uncertainties and the scatter 
into account, hardly any trend can be discerned. This is also in agreement to 
the findings for quadrupolar fluids p. 

6 CONCLUSION 

Equilibrium molecular dynamics simulation and the Green-Kubo formalism 
were used to calculate the self-diffusion coefficient, shear viscosity, and ther- 
mal conductivity for 38 different anisotropic and dipolar model fluids. A com- 
prehensive data set was obtained for each fluid and property that covers a 
substantial part of the hquid state. The statistical uncertainty of the reported 
data varies according to transport property. For self-diffusion coefficient data 
it is lower than 3 %, for shear viscosity and thermal conductivity it is around 
8 and 12 %, respectively. 

It is found that both anisotropy and dipole momentum significantly influence 
all transport properties along the bubble line. An increasing elongation always 
leads to lower values, whereas an increasing dipole momentum usually yields 
higher values. The predominant part of these variations is due to the consid- 
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erable variations in number density along the bubble line caused by L* and 
fi*"^. However, peculiar extrema are found, e.g., for the self- diffusion coefficient 
of fluids with small anisotropy. 

Density is the most important thermodynamic variable, however, compared to 
the findings for quadrupolar fluids [1], it is less dominating and the influence 
of the dipole is much stronger. 

Temperature influences all transport properties less than density. As expected, 
for higher temperatures the self-diffusion coefficient increases, the shear vis- 
cosity decreases, and for the thermal conductivity hardly any variation can be 
discerned. 



List of symbols 



D 


self- diffusion coefficient 


i 


molecule index 


j 


molecule index 


L 


molecular elongation 


m 


molecular mass 


t 


time 


T 


temperature 


u 


pair potential 


A 


increment 


At 


integration time step 


e 


Lennard- Jones energy parameter 


eo 


Electric constant 


7] 


shear viscosity 
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6 angle of nutation 

A thermal conductivity 

II molecular dipole momentum 

p density 

a Lennard- Jones size parameter 

(f) azimuthal angle 

Vector properties 

r distance vector 

Subscript 

c property at critical point 
D point dipole 

R property reduced by critical value 
Superscript 
* property reduced by molecular parameters 
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Table I 

Transport coefficients along bubble lines for 38 2CLJD fluids with different elon- 
gations L* and dipolc momentum /x*^. The numbers in parentheses denote the 
uncertainty in the last digits. 



L*=0 T* p* 


D* rj* X* 


/i*2=o 3.139 0.8062 
3.925 0.7108 
4.708 0.5838 


0.094(2) 4.65(25) 13.4(6) 
0.202(2) 2.56(16) 9.24(81) 
0.388(8) 1.48(10) 6.29(54) 


/i*2=3 3.288 0.8202 
4.110 0.7245 
4.932 0.5966 


0.083(1) 5.46(8) 13.04(89) 
0.184(4) 2.79(39) 10.12(90) 
0.372(7) 1.74(19) 7.26(73) 


/x*2=6 3.597 0.8327 
4.494 0.7345 
5.395 0.6049 


0.082(2) 6.02(21) 13.9(10) 
0.183(4) 3.14(35) 10.93(90) 
0.364(7) 1.94(3) 7.36(70) 


/i*2=9 3.954 0.8419 
4.943 0.7407 
5.931 0.6082 


0.084(2) 7.44(72) 16.41(31) 
0.187(4) 3.41(19) 13.6(10) 
0.370(7) 1.95(7) 7.01(70) 


/i*2=12 4.378 0.8456 
5.470 0.7404 
6.564 0.6024 


0.093(2) 8.99(79) 18.8(12) 
0.199(4) 3.56(20) 14.1(11) 
0.395(8) 1.62(28) 9.32(17) 


^*2=16 4.956 0.8516 
6.192 0.7400 
7.429 0.5968 


0.104(2) 18.8(16) 20.7(21) 
0.215(4) 3.75(23) 15.8(18) 
0.426(9) 1.93(7) 8.91(88) 


/i*2=20 5.505 0.8632 

6.878 0.7445 
8.254 0.5972 


0.113(2) 66.4(10) 19.9(12) 
0.226(5) 4.64(1) 14.0(20) 
0.441(9) 2.29(8) 9.93(42) 
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L*=0.2 T* p* 


D* r]* A* 


/x*2=o 2.590 0.7114 
3.237 0.6275 

3.884 0..5144 


0.088(1) 3.64(16) 11.18(62) 
0.190(1) 2.19(7) 8.7(12) 

0.300(8) 1.20(4) .5.60(35) 


/i*2=3 2.715 0.7220 
3.393 0.6367 

4.072 0.5244 


0.086(1) 4.21(42) 12.09(95) 
0.189(1) 2.19(6) 10.42(60) 
0.372(2) 1.34(1) 8.30(77) 


/x*2=6 2.952 0.7307 
3.691 0.6439 
4.428 0.5306 


0.091(1) 4.27(50) 17.7(11) 
0.197(1) 2.27(18) 15.73(86) 
0.389(1) 1.34(4) 11.3(19) 


^*2=9 3.233 0.7370 
4.040 0.6481 
4.848 0.5315 


0.098(1) 4.22(1) 26.6(16) 
0.211(1) 2.42(5) 22.0(14) 
0.415(2) 1.20(3) 11.9(19) 


11*^=12 3.572 0.7388 
4.464 0.6463 
5.355 0.5226 


0.110(1) 4.32(43) 29.4(21) 
0.235(1) 2.55(34) 24.7(31) 
0.466(1) 1.38(1) 11.5(15) 


/i*2=16 3.998 0.7438 
4.997 0.6497 
5.994 0.5229 


0.123(1) 4.75(18) 33.6(55) 
0.258(1) 2.58(23) 25.0(24) 
0.507(2) 1.18(15) 15.9(20) 
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Table I 
Continued. 



Li — U.4 1 p 


U T) A 


/x*2=0 1.899 0.5808 
2.374 0.5123 
2.848 0.4185 


0.093(2) 2.51(20) 10.0(22) 
0.187(1) 1.55(7) 7.8(11) 
0.320(3) 0.91(9) 4.84(45) 


/i*2=3 1.988 0.5881 
2.485 0.5179 


0.088(1) 2.97(15) 11.18(60) 
0.182(1) 1.80(18) 9.88(83) 

c\^'^^(o\ r\QA(<^\ <^c^^\(^f{\ 
v.oOLyzj v.ij^yo) o.yu^^ioj 


^*2=6 2.156 0.5932 
2.695 0.5225 
3.234 0.4275 


0.086(1) 2.86(17) 11.4(12) 
0.182(1) 1.84(6) 9.95(90) 
0.360(1) 0.95(2) 6.22(60) 


Ai*2=9 2.347 0.5973 
2.934 0.5254 
3.521 0.4296 


0.086(1) 3.52(18) 14.70(57) 
0.186(1) 1.86(5) 10.88(24) 
0.368(1) 0.98(3) 8.02(60) 


^*2=12 2.571 0.5988 
3.214 0.5243 
3.856 0.4257 


0.090(1) 3.49(11) 13.4(15) 
0.194(1) 1.92(19) 12.61(98) 
0.390(1) 1.16(8) 6.88(45) 
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Table I 
Continued. 



T* n KC\K T* r,* 

ij — U.OUO 1 p 


n* n* \* 


/x*2=o 1.642 0.5291 
2.052 0.4637 
2.463 0.3835 


0.092(2) 2.24(20) 9.99(84) 
0.184(7) 1.35(10) 6.01(60) 
0.303(8) 0.83(7) 4.16(35) 


^*2=3 1.724 0.5346 
2.155 0.4711 


0.087(1) 2.66(20) 9.60(79) 
0.175(1) 1.42(13) 8.10(94) 


^*2=6 1.861 0.5392 
2.327 0.4748 
2.792 0.3902 


0.082(1) 2.92(6) 10.78(16) 
0.171(1) 1.44(2) 7.57(78) 
0.334(2) 0.84(3) 5.09(45) 


Ai*2=9 2.027 0.5424 
2.535 0.4765 
3.041 0.3898 


0.080(1) 3.11(14) 10.81(69) 
0.171(1) 1.67(11) 9.7(14) 
0.342(1) 0.88(3) 6.47(99) 


^*2=12 2.207 0.5457 
2.759 0.4778 
3.311 0.3886 


0.079(1) 3.38(6) 12.58(20) 
0.175(1) 1.83(3) 10.2(16) 
0.351(2) 0.83(20) 6.92(40) 
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Table I 
Continued. 



Li — U.D 1 p 


U T) A 


/x*2=0 1.473 0.4900 
1.842 0.4321 
2.210 0.3520 


0.091(2) 2.16(30) 9.24(46) 
0.174(1) 1.14(1) 6.85(70) 
0.295(9) 0.71(2) 4.15(40) 


/i*2=3 1.536 0.4960 
1.920 0.4372 


0.082(1) 2.19(5) 7.95(31) 
0.166(1) 1.27(7) 5.41(43) 
^\'^^7^^^ c\9^'^^^o\ A7fi(A(W 


^*2=6 1.665 0.4995 

2.082 0.4398 
2.498 0.3583 


0.077(1) 2.08(39) 9.8(18) 
0.160(1) 1.34(7) 7.73(40) 
0.319(1) 0.79(10) 5.12(35) 


Ai*2=9 1.825 0.5006 
2.281 0.4378 
2.738 0.3545 


0.076(1) 2.91(12) 11.45(41) 
0.165(1) 1.75(15) 7.00(95) 
0.332(1) 0.84(1) 6.22(59) 


^*2=12 1.965 0.5062 
2.456 0.4429 
2.948 0.3603 


0.071(1) 3.28(12) 9.71(44) 
0.159(1) 1.62(4) 7.8(16) 
0.322(1) 0.89(3) 5.94(40) 
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Table I 
Continued. 



Li — U.O 1 p 


U 7] A 


//*2=o 1.230 0.4302 
1.538 0.3777 
1.845 0.3051 


0.082(2) 1.79(15) 8.28(42) 
0.157(1) 1.11(14) 4.98(18) 
0.270(5) 0.60(4) 3.86(14) 


/x*2=3 1.277 0.4364 
1.596 0.3833 
i.yio u.oiui 


0.071(1) 2.18(19) 8.52(31) 
0.145(1) 1.27(14) 6.56(50) 
^\o^i7(^\ r\f{0f7\ ^fi7(oc\\ 


/x*2=6 1.390 0.4383 
1.738 0.3835 
2.086 0.3085 


0.065(1) 2.25(45) 8.85(65) 
0.141(1) 1.27(8) 5.7(11) 
0.288(1) 0.64(4) 3.39(60) 


/x*2=9 1.513 0.4410 
1.891 0.3855 
2.270 0.3093 


0.060(1) 2.44(10) 9.1(10) 
0.137(1) 1.39(1) 7.25(70) 
0.286(2) 0.70(5) 4.32(34) 


/x*2=12 1.629 0.4466 
2.036 0.3899 
2.443 0.3155 


0.054(1) 3.02(8) 9.84(80) 
0.130(1) 1.47(2) 7.52(26) 
0.274(1) 0.74(5) 4.92(57) 
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Table I 
Continued. 



T* 1 rp* * 

Li — 1 1 p 


U 7] A 


/x*2=0 1.058 0.3970 
1.322 0.3487 
1.587 0.2860 


0.0640(5) 1.71(16) 7.12(41) 
0.128(1) 0.98(3) 4.75(35) 
0.245(1) 0.62(6) 3.50(30) 


/i*2=3 1.126 0.3976 
1.408 0.3463 


0.0590(5) 1.96(3) 7.33(39) 
0.127(1) 1.18(6) 4.73(13) 

^\of{7(^\ c\'^f{(A\ Ar\Q(^Q\ 
v.zotyij v.ooy^) ^.yjijy^oiyj 


^*2=6 1.207 0.4041 
1.508 0.3524 
1.810 0.2887 


0.0480(5) 2.50(9) 7.45(78) 
0.113(1) 1.24(20) 6.04(29) 
0.232(1) 0.67(4) 3.62(17) 


Ai*2=9 1.308 0.4101 
1.635 0.3576 
1.961 0.2900 


0.0390(5) 3.31(13) 8.54(80) 
0.104(1) 1.56(15) 5.22(18) 
0.228(1) 0.72(5) 4.30(25) 


^*2=12 1.423 0.4163 
1.779 0.3625 
2.134 0.2923 


0.0330(5) 4.09(10) 8.98(30) 
0.0970(5) 1.75(5) 6.64(98) 
0.223(1) 0.75(1) 3.92(34) 
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List of Figures 



1 Phase diagrams for two selected elongated 2CLJD fluids 
(L*=0.2) where one is non-polar and the other strongly 
dipolar. Saturated densities, taken from [20], are represented 
by the lines joining at the critical point depicted by •. The 
investigated state points are indicated by o for /i*^=0 and by 
A for /x*2=12. 

2 Saturated liquid density of Stockmayer fluids (L*=0) 
for three reduced temperatures (Tr=0.6, 0.9, and 1) as 
functions of the dipole momentum. The solid lines indicate 
the overall correlation in [8], whereas the symbols are 
from correlations fltted to each model fluid's vapor-liquid 
equilibrium individually. 

3 Self-diffusion coefficient of spherical (L*=0, empty symbols), 
elongated (L*=0.505, grey symbols), and strongly elongated 
(L*=l, full symbols) 2CLJD fluids over reduced density along 
bubble lines. Reduced temperatures vary from Tr=0.6 to 0.9. 
Lines are guides for the eye. 

4 Self-diffusion coefficient of spherical (L*=0, empty symbols), 
elongated (L*=0.505, grey symbols), and strongly elongated 
(L*=l, full symbols) 2CLJD fluids over number density along 
bubble lines. Reduced temperatures vary from Tr=0.6 to 0.9. 
Lines are guides for the eye. 
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Self-diffusion coefficient of spherical (L*=0, empty symbols), 
elongated (L*=0.505, grey symbols), and strongly elongated 
(L*=l, full symbols) 2CLJD fluids over number density in the 
homogeneous liquid at rR=0.9. Lines are guides for the eye. 

Self-diffusion coefficient of spherical (L*=0, empty symbols) 
and strongly elongated (L* = l, full symbols) 2CLJD fluids over 
reduced temperature in the homogeneous liquid along different 
isochores. o: p*=0.8062. A: p*=0.8327, 0: p*=0.8456. Lines are 
guides for the eye. 

Shear viscosity of spherical (L*=0, empty symbols), elongated 
(L*=0.505, grey symbols), and strongly elongated (L*=l, full 
symbols) 2CLJD fluids over number density along bubble 
lines. Reduced temperatures vary from Tr=0.6 to 0.9. Lines 
are guides for the eye. 

Shear viscosity of spherical (L*=0, empty symbols), elongated 
(L*=0.505, grey symbols), and strongly elongated (L*=l, 
full symbols) 2CLJD fluids over number density in the 
homogeneous liquid at Tr=0.9. Lines are guides for the eye. 

Shear viscosity of spherical (L*=0, empty symbols) and 
strongly elongated (L*=l, full symbols) 2CLJD fluids over 
reduced temperature in the homogeneous liquid along different 
isochores. o: p*=0.8062. A: p*=0.8327, 0: p*=0.8456. Lines are 
guides for the eye. 
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10 Thermal conductivity of spherical (L*=0, empty symbols), 
elongated (L*=0.505, grey symbols), and strongly elongated 
(L*=l, full symbols) 2CLJD fluids over number density along 
bubble lines. Reduced temperatures vary from Tr=0.6 to 0.9. 
Lines are guides for the eye. |36 

11 Thermal conductivity of spherical (L*=0, empty symbols), 
elongated (L*=0.505, grey symbols), and strongly elongated 
(L*=l, full symbols) 2CLJD fluids over number density in the 
homogeneous liquid at Tr=0.9. Lines are guides for the eye. |37 

12 Thermal conductivity of spherical (L*=0, empty symbols) and 
strongly elongated (L*=l, full symbols) 2CLJD fluids over 
reduced temperature in the homogeneous liquid along different 
isochores. o: p*=0.8062. A: p*=0.8327, 0: p*=0.8456. Lines are 
guides for the eye. l38 
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Fig. 1. 
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Fig. 2. 
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Fig. 3. 
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Fig. 4. 
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Fig. 5. 
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Fig. 6. 
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Fig. 7. 
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Fig. 8. 
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Fig. 9. 
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Fig. 10. 
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Fig. 11. 
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Fig. 12. 
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